The s— wave implosion and explosion of free particles 
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A free quantum particle in two dimensions with vanishing angular momentum (s-wave) in the form 
of a ring-shaped wave packet feels an attractive force towards the center of the ring, leading first to 
an implosion followed by an explosion. 
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The spreading of a wave packet is an important topic 
in quantum mechanics It arises due to different mo- 
menta contained in the packet. As pointed out by Max 
Born [| even a classical ensemble of particles spreads. 
The only quantum effect lies in the connection between 
the position and the momentum distributions, that is, 
the spreads in position and momentum are not indepen- 
dent. In the present paper we predict a wave effect of 
the free particle that manifests itself in the dynamics: A 
ring-shaped wave packet in two dimensions of vanishing 
angular momentum feels an attractive force towards the 
center of the ring. This effect is unique to two dimensions 
and manifests itself most strongly for vanishing angular 
momentum. 

In order to illustrate this surprising phenomenon we 
first consider the free motion of wave packets in N di- 
mensions. We then concentrate on two and three dimen- 
sions, where we propagate a circular shell and a spher- 
ical shell, respectively. All wave packets have initially 
vanishing angular momentum and vanishing average ra- 
dial momentum. We show that the circular shell first 
implodes and then explodes. In contrast, the spherical 
shell does not implode but only explodes. We demon- 
strate this wave effect unique to two-dimensional space 
using (z) analytical expressions for the average radial po- 
sition and momentum, (ii) the Wigner function approach 
1^ towards quantum mechanics, {Hi) analytical expres- 
sions Q for the radial Green's function in two dimensions 
and (iv) the numerical propagation of wave packets. 

Our analytical approach considers the propagation of 
the initial wave function 
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in the absence of any potential. Here M = 
[2/F(iV/2-|-7)]^/^(5r-^/2-'r denotes the radial normal- 



ization constant and r 



) is the ra- 
dial variable in TV space dimensions. Moreover, 5'^^-' = 
27r^/2/p (•^/2) is the total soHd angle. The quantity 6r 
defines a characteristic length scale of the packet and is 
a measure of the radial width H . 



We first take the exponent 7 to be equal to 2. This 
special form of the initial wave function - polynomial 
in coordinates multiplied by a Gaussian - allows us to 
calculate analytically the time dependence of the average 
radial position (r^^^ ) and the average radial momentum 
(p^^-') in TV dimensions. We then turn to the general case 
of arbitrary 7 > 1 but concentrate on two dimensions. 

In the derivation of our results we follow three inde- 
pendent approaches that provide identical results: (i) We 
consider the time evolution in Wigner phase space and 
calculate the averages with the help of the Wigner func- 
tion; (ii) we use the method of the generating function 
Q to propagate a modulated Gaussian wave packet of 
the free particle in N dimensions, and evaluate the av- 
erages; {Hi) we propagate the free wave packet with the 
help of the radial Green's function of the free particle 
in two dimensions. Here we do not present the detailed 
calculations but only summarize the results 

We start by presenting first the general results for 
the propagation in N dimensions and then focus on 
the special case of two and three dimensions. It is 
convenient to introduce the dimensionless time variable 
T = ht/{Sr'^M), where M denotes the mass of the parti- 
cle. We find the expression 
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for the average radial position. Here we have introduced 
the abbreviation 



x(^) = 1 + 



and the initial radial position 



(3) 



(4) 



in N dimensions, where F denotes the Euler gamma func- 
tion. 

Equation (||) shows that for large times the average ra- 
dial position increases linearly with time, corresponding 
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to an explosion of the wave packet. This holds true for 
any number of dimensions. 

For short times we can expand the square root in Eq. 
(0) which yields 
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An interesting situation occurs when a'^^ > 2. In this 
case, the average radius initially decreases, corresponding 
to an implosion of the wave packet. From the definition 
Eq. (H) of a^^-* we arrive at the implosion condition 



(iV-l)(iV-3) <0. 
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Hence an implosion occurs only for N ~ 2 corresponding 
to a(2) = 15/7 = 2 + 1/7. 

The phenomenon of the implosion can also be viewed 
by considering the time evolution of the average radial 
momentum 

(p<^^>(r)-M|(rW)-i(r(^))i. (7) 
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determined by the time derivative of the average position 
(r'^^^), which yields 



where 
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The radial momentum is negative, corresponding to 

1 /2 

an implosion, provided r < [a^^^ — 2] . Since in two 
dimensions a^^^ — 15/7 this condition reads r < r^in = 
1 / \/7. At that moment the momentum vanishes and the 
implosion turns into an explosion. The average radial 
position assumes a minimum rmin = (f^^^) (Tmin)- 

In order to highlight the special property of two dimen- 
sions we now compare the dynamics of the wave packet 
Eq. (|l|) for 7 = 2 in two and three dimensions. In Fig. |l| 
we show the time evolution of the average radial position 
and the average radial momentum in two (solid curves) 
and three (dashed curves) dimensions based on Eqs. (H) 
and d) with a(2) = 15/7 and a^^) = 2. 

For short times the average radius of a ring-shaped 
wave packet decreases as a function of time. This initial 
implosion of the circular wave is followed by a final ex- 
plosion. We emphasize that this phenomenon is not due 
to an external classical force but occurs for a free particle 
moving in two dimensions Obviously the spreading 
of the ring-shaped wave packet in the radial direction 
is asymmetric; initially the probability flux towards the 
center is larger than towards the outside. 

In the three-dimensional case, where we propagate a 
spherical shell, the average position and the average mo- 
mentum do not display the soft implosion effect of two 



dimensions. Indeed, since a^^^ = 2, the average momen- 
tum is always positive and the average position cannot 
decrease but must increase. 

The contraction effect is maximal at Tmin = l/v^ and 
according to Eq. (|^) the average radial position 
assumes the minimal value 
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Hence, it is a rather soft implosion. 

We note that this result is independent of 5r, since 
there is only one length scale in the problem. When we 
choose a different form for the packet we can make the 
effect larger. For example, by replacing in Eq. (|^) by 
sin(r^/(Sr^), we obtain at Tmin — 1-11 
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In order to study the implosion effect we concentrate 
on two dimensions and now consider an arbitrary value 
of the exponent 7 > 1 in the the wave packet Eq. 
With the help of the corresponding radial propagator 
we obtain the radial wave function which yields the ex- 
pression 
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for the average radial position. Here the integral 
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containing the confluent hypergeometric function iFi can 

be expressed in terms of the generalized hypergeometric 

function 3^2. However, we found it more convenient to 

evaluate this integral numerically. For each value of 7 we 

f 21 

obtain an implosion but the ratio rmin/rQ is always of 
the same order of magnitude. 

We now turn to the explanation of the implosion effect 
using the Wigner function ||^. The free particle is one of 
the few quantum mechanical systems in which the clas- 
sical phase space dynamics is identical to the quantum 
mechanical one. Indeed, the Wigner function satisfies the 
classical Liouville equation. Hence, we find the Wigner 
function at time t from the Wigner function at time t = 
by replacing the position fhy r — pt/AI. We emphasize 
that now r and p denote A'^-dimensional vectors. The 
corresponding probability distributions in position and 
momentum follow by integration over the conjugate vari- 
able. 

Since for the free particle classical and quantum phase 
space dynamics are identical the quantum effect of im- 
plosion rather than explosion is stored in the initial 
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Wigner function. To be more precise, it is in the cor- 
relations between position and momentum. A classical 
phase space distribution enjoys a probability interpreta- 
tion and therefore always has to be positive. In contrast, 
the Wigner function can assume negative values. Indeed, 
both the Wigner functions of the circular shell and of the 
spherical shell contain domains in phase space where they 
become negative. Since the Wigner function is normal- 
ized the total volume of the Wigner function is unity. 

(2) (3) 

We can calculate numerically the volumes V_ and V_ 
corresponding to the negative parts. We find in the two 

(2) 

dimensional case V_ ~ 0.27, whereas in three dimension 

VI ' ~ 0.23. Hence, there is a slightly larger contribu- 
tion of negative parts in two dimensions than in three 
dimensions. 

In the phase space description of quantum mechan- 
ics quantum effects are stored in the negative parts of 
the initial Wigner function. Since a classical ensemble of 
particles cannot have a negative phase space distribution 
it cannot display the implosion effect. Needless to say 
a classical ensemble with a ring-shaped initial position 
distribution can show an implosion provided there exists 
a non-vanishing negative initial radial momentum. How- 
ever, in our quantum mechanical wave packet the average 
radial momentum vanishes. 

For the choice Eq. (|l|) of the initial wave packet the 
radial probability density 

= I |2 r^+27-ie-(rMr)^rf^ (14) 

depends on the number of dimensions. Hence, one might 
argue that the implosion is a result of the slightly differ- 
ent initial radial wave functions. In order to exclude this 
argument we have performed a numerical integration of 
the Schrodinger equation 
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for the radial wave function Wg^'' of the free particle in 
two and three space dimensions, corresponding to zero 
angular momentum. Here we always start with the same 
initial radial wave function. In this way we guarantee 
identical radial probabilities densities in two and three 
dimensions. 

The radial potential 



u^^\r,t) (15) 
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in N dimensions involves the product {N — 1){N 
familiar from the implosion condition Eq. (^). 

In three dimensions the potential V^'^'' vanishes. How- 
ever, for N = 2 the quantum anti-centrifugal potential 



(17) 
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does not vanish and is negative, that is, attractive 
This potential is responsible for the implosion. 

A large class of wave packets with vanishing angu- 
lar momentum displays this effect of first imploding and 
then exploding ||l^. Indeed, in the insert of Fig. |l] we 
show the short time evolution of the average radial posi- 
tion and the average radial momentum in two and three 
dimensions for an initial radial Gaussian wave packet 



displaced by an 



u{r,t = 0) = A/'cxp — {r ~ p) /5r'^ 
amount p from the origin. As in the analytical results 
we find an implosion of a circular but not of a spherical 
shell. 

So far we have concentrated on a single particle but 
the effect of implosion manifests itself also in the relative 
motion of two free particles. We can envision the total 
wave function of the two particles to be of the form Eq. (|^) 
corresponding to an entangled state. When we assume 
that each particle is allowed to move in a D-dimensional 
space the total wave function lives in a = 2D dimen- 
sional configuration space. The physical situation of a 
single particle imploding in two dimensions corresponds 
to two particles constrained to one dimension. When 
they move along a line their average relative separation 
first decreases and then increases. It looks like the par- 
ticles first attract and then repel each other. However, 
when they move in two or three dimensions, their relative 
distance always increases. 

We conclude by briefly outlining an experiment to 
demonstrate the implosion effect. At first sight, one 
might think of cold atoms around a wire familiar from 
atom chips [ p^ . However, it might be easier to make use 
of the analogy between the Schrodinger equation and the 
paraxial wave equation of classical electrodynamics. We 
shine monochromatic light onto an opaque screen with 
an annular aperture. In order to avoid diffraction from 
sharp edges, we introduce apodization, that is, a posi- 
tion dependence transmission coefficient. At the edges of 
the ring almost no light is transmitted, whereas at the 
center of the ring almost all light passes through. This 
arrangement creates an intensity distribution after the 
screen confined mainly to a ring. The subsequent free 
propagation in space is analogous to the free propaga- 
tion in time of the Schrodinger wave function. The time 
dependent probability distribution then corresponds to 
the transverse intensity distributions at different loca- 
tions downstream from the screen [ p^ . 
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S. is partially supported by DFG. ready when we approximate the radial propagator in 

two dimensions by the difference of two one-dimensional 
Green's functions of free motion. Hence, the implosion 
effect is a modification on top of the Poisson spot. 
[14] M. Born and E. Wolf, Principles of Optics (Cambridge 
University Press, Cambridge, 1999) p. 417. 
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FIG. 1. Comparison between the time evolutions of a cir- 
cular (solid line) and a spherical shell (dashed line) using the 
average radial position (top) and the average radial momen- 
tum (bottom) obtained from Eqs.(|^ and ^ for a'^' — 15/7 
and a'^*) = 2. On average the two-dimensional shell first con- 
tracts and then expands, whereas the three-dimensional one 
expands from the onset. In the left inserts we show the short 
time behavior of a displaced Gaussian radial wave packet ob- 
tained numerically. Here we have used the same Gaussian for 
an initial radial wave function in two and three dimensions 
with Sr = 0.4 and p = 1.5 in arbitrary units. Position is 
scaled in units of the initial position Tq'^' of the wave packet, 
and momentum in terms of the final momentum . 
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